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Abstract 

We study the CP-violating, |A5| = 1 parameter e'/e by computing the hadronic 
matrix elements in the chiral quark model. We estimate in the chiral expansion the 
coefficients of the next-to-leading order terms that correspond to the operators 
and Qii- We consider the impact of these corrections on the value of e'/e. We also 
investigate the possibility that the chiral quark model might drastically reduce the 
dimensional regularization 75-scheme dependence of current evaluations of e'/e. 
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1 Foreword and Summary 

A detailed study of the CP-violating, \AS\ = 1 parameter e'/e puts us at the 
crossroad of different aspects of standard model physics and, accordingly, provides a 
promising testing ground of model-dependent features. 

As the precision in the experimental measurements of e'/e improves all theoret- 
ical predictions face the challenge of reducing their uncertainties. These have a twofold 
origin as, on the one hand, the analysis of the Wilson coefficients must be pushed to 
higher orders in as and a (as well as including all the relevant operators) while, on 
the other hand, the estimate of the hadronic matrix elements has to be substantially 
improved. 

Indeed, a great deal of work has been done in this direction. In recent years the 
Wilson coefficients of the effective lagrangian for \AS\ = 1 weak decays have been 
computed to the next-to-leading order (NLO) ||2|, |3|, and matrix elements have been 
calculated by a variety of techniques like 1/Nc @, quark models |^] and lattice Yet 
it is perhaps fair to say that, all efforts notwithstanding, the current uncertainty in the 
theoretical prediction of e'/e is still large, most of it arising from the hadronic matrix 
elements. 

There are at least four relevant perturbative expansions which are important in 
estimating hadronic matrix elements. These are the 1/Nc and a^/vr expansions for 
non-factorizable contributions, which have been studied in [Q] and Q, respectively. 
Then we have the chiral expansion, which comes in two varieties according to the 
energy scale at which it is considered. The coefficients of the chiral lagrangian take 
contributions from chiral loops (corresponding to momenta below the scale Aqcd at 
which all quark fields are integrated out) as well as from light-quark loops (between 
Aqcd and the chiral symmetry breaking scale A^). 

The chiral quark model (xQM) |^, 0, ^ can be thought as interpolating between 
the lower-end of the perturbative short-distance analysis (~ 1 GeV) and the chiral 
lagrangian energy region (Aqcd ~ M). In the model, the scale M is identified with the 
constituent quark mass. This is the only free parameter in our analysis and most matrix 
elements are a sensitive function of it. What are the allowed values for this parameter? 
As we shall discuss in section 2, a rather large range (100 < M < 400 MeV) is possible, 
different analyses having different outcomes according to the actual physical processes 
considered. 

In this paper, we begin a study of e'/e in the xQM by focusing on two aspects: the 
relevance of the NLO corrections to the hadronic matrix elements and the problem of 
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the 75-scheme independence of the estimate. Our original motivation was to study the 
gluonic magnetic dipole penguin operator 

Y£2m,S(T^,t'^Gr(l-75)'i, (1.1) 

which until recently has been tacitly neglected in all discussions of e' je. The potential 
relevance of the operator Q\\ for e' je has been discussed in a previous work g. We 
confirm the result of ref. |^ up to an overall sign. In section 3 we discuss in detail 
the bosonization of Q\\. The first non- vanishing contribution of Q\\ to the K — > vrvr 
matrix elements comes from a term that is NLO in the chiral expansion (©(p"^)). Even 
though the corresponding coefficient of the chiral expansion — as well as the short- 
distance Wilson coefficient — is large, the matrix element turns out to be kinematically 
suppressed by a factor m^/m|^ with respect to the dominant operators, making its 
contribution to e' je sub leading. 

In order to better gauge the significance of this NLO contribution, we estimate, in 
section 4, the NLO correction to the matrix element of the ordinary gluonic penguin 
operator 

which gives the dominant contribution to e' je. Such an NLO correction to the hadronic 
matrix element of Q% turns out to be generally larger than the contribution of the 
operator Qn, as well as the leading-order (LO) contributions of the operators Qs, Q5 
and Q-j (see section 5 for definitions). Its size increases monotonically with M, ranging 
from 14% of the LO term for M = 120 MeV up to 25% and more for M 220 MeV. 
Our perturbative results (we perform an expansion in powers of MVA^) should not be 
trusted for larger values of M, higher-order corrections being needed. 

Our approach makes it also possible to study the dimensional regularization 75- 
scheme dependence of the matrix elements of these operators. This is an important 
problem since it causes, if we take the NLO results for the Wilson coefficients in the two 
schemes, a significant part of the current theoretical uncertainty (see for instance ref. 
H) — as large as 70-80% for rat = 170 GeV. One expects this scheme dependence to 



disappear when matrix elements are consistently evaluated [10|. This cannot be done 
within the lattice or 1/A^c approach. 

As a preliminary analysis, in section 5 we use the xQM to study the 75-scheme 
dependence of e' je in a toy model, made of the LO matrix elements of the operators 
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Qq (for the / = amplitude) and of the electroweak penguin operator (for 1 = 2): 
3 

Qs = 2 S«7m(1 - 75)d/3 ^1 9/37''(l + J5)qa , (1-3) 

Q 

which becomes important for a heavy top quark. As a matter of fact, the Wilson coef- 
ficients of these two operators exhibit the largest scheme dependence. In comparison, 
the scheme dependence of the Wilson coefficients of (^4 and Qg, next in relevance to 
Qq and Qs, is ten times smaller. 

We find that the 75-scheme dependence of e' /e in our toy model is reduced below 
the 10% level for values of the parameter M well within the allowed range. This means 
that, for any value of ^qcd in the 200-400 MeV range and of the renormalization scale 
H between 1 GeV and rric, it is possible to find a value of M in a restricted range (120- 
220 MeV) for which the scheme dependence of the hadronic matrix elements exactly 
cancels that of the Wilson coefficients. 

Notice that in the xQM the matrix element of Qn turns out to be scheme-indepen- 
dent, as it happens for its "LO" Wilson coefficient |11]. 

The range of values for which we find scheme independence is consistent with the 
values of M for which we can trust our perturbative results at the order M^/A^ 
(M ^ 200 MeV). To extend our analysis in a reliable manner to values of M of 300 MeV 
and more requires considering higher-loop effects — a difficult computational problem. 

As encouraging as these results might be, we should remember that we are con- 
sidering a toy model for e' /e in which only the two leading operators Qq and Qg are 
included. In order to ascertain conclusively the scheme independence and the size of 
e'/e, all relevant operators should be evaluated in the xQM and consistently included. 
Preliminary estimates show that such an extension stabilizes the range of values of 
M for which the scheme independence is substantially weakened, further reducing the 
renormalization scale dependence. 

In the final part of the paper (section 6) we compare the contribution to e' /e of the 
dipole penguin operator Qu to the other ten traditional operators. Since the hadronic 
matrix elements of the other operators beside Qq are known in the only to the 

leading factorizable order, we follow for the ten traditional operators the 1/A^c approach 
ofref. i. 

A detailed anatomy of all individual contributions is summarized in tables 6, 7 and 
8. We find these tables a useful way of displaying the available information and a 
standard reference for future developments. 

The Wilson coefficients are computed to the NLO level, at least to the extent made 
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possible by the entries of the anomalous dimensions currently available — the three- 
loop NLO mixing between Qu and the other operators and related two-loop initial 
conditions are still missing. This may amount to a 40-50% uncertainty in the Qu 
Wilson coefficient at the 1 GeV scale (extrapolating the results for 6 — > 57 given in 
ref. [1l2| ). Nonetheless, this little affects the final result for e' /e since the evolution of 
the first ten Wilson coefficients is independent on the presence of the Qu (and Q12) 
operator. 

Our analysis shows that the inclusion of the magnetic penguin operator affects the 
prediction of e'/e always below the 10~^ level, therefore putting Qu in the same class 
as other operators like Q3, Q^ and Qi, whose contributions are even smaller; yet, for 
mt ^ 170 GeV, the size of e'/e itself becomes of the same order as the Qu contribution. 

The completion of the study of e'/e in the xQM approach is needed to ascertain 
conclusively the 75-scheme independence of the result and, what is equally relevant, the 
corresponding size of e'/e. Corrections of the order of the current quark masses and 
a/vr will then have to be included. We expect the mass corrections to be of the order 
or less than 10%, while those of (soft) gluon contributions to be larger ( ^ 50%). We 
recall that in our toy model we find systematically a factor of 3-4 enhancement with 
respect to the average of the 1/Nc results with Qq and Qs only and mt = 170 GeV. 
This is partially due to a reduced cancellation between the two operators in the xQM. 

We also stress the relevance of the calculation of NLO corrections to all relevant 
matrix elements, which become crucial for mt ^ 170 GeV, due to the large cancellations 
that occur between gluonic and electromagnetic penguin contributions. 



2 BOSONIZATION OF QUARK OPERATORS 

In general, the effective lagrangian at the quark level has the form 
£,Q = J2CiQ,, (2.1) 

i 

where the coefficients Cj contain the short-distance physics above the renormalization 
point ^ ~ 1 GeV, and Qi are operators containing quark fields (at the scale fi). When 
the quark operators are bosonized, they will be represented by some linear combination 
of meson operators jCj^^: 

Qi^Y. ' (2-2) 

j 
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leading to the following chiral lagrangian at the meson level: 
where Gj{Qi) = Ci bij (no sum on i). 



One way to perform the bosonization is to give the complete list of oper ators C'f'^ to 



which some quark operator Qi could a-priori contribute |13, 14, 15]. To determine the 
various coefficients bij, one tries combinations of experimental data and phenomeno- 
logical estimates. In this work, we have chosen a version of the xQM advocated by 
various authors |^, |^ . This type of effective low-energy model is obtained by adding a 
new term to the ordinary QCD lagrangian: 

C^QM = -M(qiiJ:qL + qL^kR) , (2.4) 

where q = (u,d,s), and the 3x3 matrix S = exp(2in//) contains the pseudoscalar 
octet mesons 11 = X^a'^'^^"/^ (a = The scale / is identified with the (bare) 

pion decay constant, so that numerically / — /tt — 93.3 MeV (and equal to fx in the 
chiral limit). 

The parameter M is characteristic of the xQ,^ ^-nd represents a typical constituent 
quark mass for the light quarks. For M ^ the chiral symmetry is broken and the 
quark condensate acquires a non-vanishing vacuum expectation value. A central value 
of 320 MeV is found in ref . |p , albeit with a comparable uncertainty. Smaller values are 
found by fitting data relevant to meson physics, as in ref. |l^, where M = 223 it 9 MeV 
is quoted. It seems fair to consider the range 100 MeV < M < 400 MeV, as a reasonable 
present-day uncertainty on the determination of the parameter M. 

The resulting field theory with pseudoscalar quark-meson couplings (see appen- 
dix A) and quark loops gives a good description of the K ^ -kit amplitudes to leading 
order in 1/A^c and under certain conditions of the anomalous vr'' — > 77 amplitude 
1 17, |l^, as well as — > 77 [18, 19|. It is also worth recalling that the A = 1/2 rule 
in CP-conserving kaon decays is reproduced within the model up to factors 2-3 [Q, 
whereas other approaches fail by an order of magnitude. 

When performing loop integrals within this model, the logarithmic and quadratic 
divergences are identified with /^r and the quark condensate {qq), respectively, that is: 

/i°^ = l^/- (g5)(°) = ^C,, (2.5) 

where and {qq)^^^ are the pion decay constant and the quark condensate to ze- 
roth order in ag (no gluon-condensate corrections and zero current-quark masses, see 
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eq. ( 4.20D ). We distinguish between fi^^ and / because of their different origin, even 
though at out level of approximation fi^^ = f — fw 
In dimensional regularization (d = 4 — 2e) we have 

,2 



= ^-7i? + ln47r + ln^ 



and 



Cg = M'if^ + 1) 



(2.6) 



It is perhaps useful to recall that, introducing in the effective theory a cut-off A one 
obtains [Si 



U = log(AVM2) + 



and 



Cq = -A2 + log(AVM2) + 



(2.7) 



To the order we work we cannot remove the cut-off. Equations ( |2.5| ) are thus a conve- 
nient book-keeping device to identify f-j^ and (qq) in the loop integrals. 

The model has a "rotated" picture, where the term C^qm in (2^) is transformed 
into a pure mass term —MQQ for rotated "constituent quark" fields Ql r '■ 



QL 



Ql = iqL and q^ ^ = ^^q^ , 



(2. 



where ^ • ^ = S. The meson-quark couplings in this rotated picture arise from the 
kinetic (Dirac) part of the constituent quark lagrangian. These interactions can be 
described in terms of vector and axial vector fields coupled to constituent quark fields 
Q = Qr + Ql-- 



where 



Q 



7^V^ + 7^75A 



Q: 



(2.C 



and 



(2.10) 



(2.11) 



Here and are the external fields containing the photon (as well as the W field). 

Using (p.Sj), the strong chiral lagrangian 0{p'^) can be understood as two axial 
currents attached to a quark loop, leading to 

4^) ~ Tr U^^H . (2.12) 



Using the relations |16, 20 1 



2iA, 



(2.13) 



6 



one obtains the leading strong chiral lagrangian 



f2 

42) ^ i_ Tr (Z)^StL>^S) , 



(2.14) 



where Z)^ is the covariant derivative. Note that is invariant under local chiral 
transformations in agreement with the invariance of £s • In contrast, the vector 
field transforms as a gauge field. Therefore, terms involving the vector field would 
manifestly break local chiral invariance and cannot appear in a chiral lagrangian, except 
for the anomaly term [|l9|. The vector field can only be present indirectly through 
some field tensor obtained from the commutator of the total covariant derivative of 
the lagrangian. In addition, the covariant derivative of may occur in higher-order 



terms P, 15, 2C, 21|. 



In the version of the xQM given by ( |2.4[ ) (the unrotated picture) the momenta which 
correspond to derivatives of the fields have to be extracted from the quark propagators 
in the loop diagram, and it is in general difficult to see the correspondence between the 
quark and the bosonized version of the operators. The rotated picture, where the axial 
vector fields in eq. (^) couple to the quark loops with derivative couplings, is more 
transparent in this sense. For example, for 0(p'^) terms (with two momenta and one 
current quark mass), no momenta have to be extracted from the quark propagators, 
and we can readily deduce the chiral lagrangian (bosonized quark operator) in terms 
of the rotated fields. 



3 BOSONIZATION OF Qn 



The magnetic dipole operator ( |1.1D can be rewritten as 
9s 



Q 



11 



87r2' 



mdR + nis L 



a ■ G d + h.c. 



(3.1) 



where a- G = a^^ Gl^^t"", G'^,^ is the gluon field, are the SU(3) generators, normalized 
as Tr{tH'') = (l/2)(5"^ and {L,R) = (1 =F 75)/2. The Hermitian conjugate is there to 
remind us of the corresponding s — > d transition. This will be understood throughout 
the paper. The dipole operator in eq. (3.1) can be written in a chiral SU{3) covariant 
form as 



Qii = A Ur -^9^- <^-GqL + qL (^-G X-Ml 



8tt 



QR 



(3.2) 



where M.^ 



diag (m„, m^;, m^) is the current quark matrix. The Gell-Mann ma- 



trices X± = {Xq ± iX-j)/2 project AS* = ±1 transitions out of the quark fields 
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Qi = {ul, cIljSl)- Note that in eq. ( |3.2D the operator transforms as {8i,lji) un- 
der the chiral SU{3)l x SU{3)r symmetry if the current quark matrix is taken to 
transform as Aiq UiiAiqUj^, where Uji and Ul are the chiral SU{3) transformation 
matrices. 

In the rotated picture, Qu reads 



Q 



11 



8^2 



where F, 



V.A 

(-) 



cT-GQ, 
/2, with 



and 



(3.3) 



(3.4) 



The operator corresponding to Qu is understood in terms of a quark loop. Let us 
first note that the lowest-order contribution is given by a quark loop where only the 
field F^-^ is interacting 



£^^\Qu) ~ Tr(F(^)) = Tv\^^MqX- + X-Ml^ 



(3.5) 



This is the contribution that cancels against the pole term induced by a non-vanishing 
{0\Qii\K^) transition [22|, in agreement with the FKW theorem [23|. To the NLO, we 
obtain a term corresponding to an interaction of F^--^ and two axial fields attached to 
a quark loop: 



C^'^HQu) ~ Tr 



(3.6) 



Note that F^^ is not contributing in eqs. (|3.5| )-( |3l6[ ), because there must be an even 
number of 75 's in the quark loop. Using (|2.13|) and (0), we find that £^^\Qii) can 
be written in the form 

C^^HQii) = (Qii) Tr [(stAl,A_ + A^AlJs) Df^J^^D^T.] . (3.7) 

In eq. ( |3.2D , we see that A_ and Mq are always next to each other, and so they must 
be in the bosonized version of Qu- We neglect the possible terms appearing as the 
product of two traces. The latter correspond to two quark loops connected by (soft) 
gluons and will therefore be suppressed by a factor (as/vr)^, at least. Due to the CPS 
symmetry [24| of the quark operators, also the corresponding bosonized operators have 
to exhibit such a symmetry. As a consequence the two terms in eq. ( |3.7D must appear 
with equal weight. Among all the possible terms for \AS\ = 1 kaon decays, which can 
be obtained by inserting the flavor-changing matrix A_ into various places within the 
trace in ci^^ [15|, we are thus left with the unique term in eq. ( |3.7| ). 
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The coefficient in eq. ( p.TD can be most easily calculated in the unrotated picture, 
by considering, for instance, the off-shell X — > tt transition. In order to calculate 
G'q'\Qii), we expand to find the amplitude 



A{K+ ^7T+;Qu) = {7r+\i C^iQu)\K+) = ^ (m, + ma) (Qu) , (3.8) 



where k is the off-shell momentum. Matching eq. ( |3.8[ ) with the corresponding quark 
loop amplitude, represented by the diagrams in Fig. 1, leads to 

Gt\Qn) = -'4im)G(3\) , (3.9) 



4 '''' Vievr . 
where we have used 

which represents the two-gluon condensate contribution to the quark condensate (sec- 
ond diagram in Fig. 2). Diagrams analogous to those in Fig. 1 with the two meson 
lines entering at the same point do not contribute to momentum dependent terms, 
and are disregarded. In Table 1 we show the values of {qq)G as a function of M. 
For the gluon condensate we have taken the central value of the lattice evaluation 
i^GG) = (460 ± 20)^ MeV^ The entries shown can be scaled accordingly for 



Having determined G^\Qii) from the — > tt transition, we can deduce the K 



other values of the gluon condensate. 

Having determined G^\Qii) fron 
2tt amplitude from the chiral structure of C^^^Qu): 

^(i^o^7r+7r-;Qii) = ^ (m, - m^) (Qn) , (3.11) 



Even though the coefficient in C^^\Qii) is large, the Qu contribution to X ^ 27r 
is small because of the factor (in place of obtained for Qq). The modest role 
played by Qu in e'/e comes therefore from this kinematical suppression rather than 
from its contribution being NLO in the chiral expansion. 

4 Next-to-leading order bosonization of Qq 

The bosonization of Qq follows basically the same line as for Qu in the preceding 



section. The standard expression (obtained from eq. (1.2) by a Fierz transformation) 



Qe = -8Y.isLqR)iqRdL) , (4.1) 
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can be rewritten in the rotated picture as 

Qe = - 8 (QL)a(.QR)5(QR)s{QL)f3 , (4.2) 

where = ^A_^^, and the greek letters are flavor indices. Thus, hkewise to 

£^^\Qii), the chiral representation of Qq to leading order can be written as 



Tr 



(4.3) 



which by means of eq. ( p^ ) can be written in the same familiar form as the other 
IAS" = 1| octet operators 0{p^) |2|]: 



(2) 

AS=1 



(Qe) Tr(A_Z)^StD^E), 



(4.4) 



where G^\Qq) is the most important component. The term in eq. (4.4) gives rise to 
the amphtudes 



72 ^'4 (^6) 



^(i^o^vr+vr-; Qg 



V2 

P 



2 2 



(4.5) 
(4.6) 



Now, we want to find the bosonization of Qq to the same order as C^'^\Qii) in the 
preceding section. That is, the bosonized operator has to contain a mass insertion in 
addition to what is already included in C'^'^\Qq). The QCD mass lagrangian can be 
written as 



C 



which can be transformed to the form C 



QMqQ, where 



Therefore, a possible NLO representation of Qq is given by 
£(4)(Q6) ~ Tr [F(_)^^Mf 



(4.7) 



(4.^ 



(4.9) 



with the addition of two other terms, where the quantities within the trace are per- 
muted. In eq. ( ^ ) represents some part of Mg (i.e. X = R, L,V, A). Using 
eq. ( ^.13| ), we obtain three different CP5'-symmetric terms: 



£g)(Q6) = Gg^(Q6) Tr [(stA^,A_ + X.MI^) Df^^^D^^ 
^%\Q^) = G%\Qq) Tr [(A_StAlg + ^ A_) D^'T^D^ll 



(4.10) 
(4.11) 
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4^(Q6) = G^^(Q6) TV 



(4)/ 



(4.12) 



These three terms are in fact hnear combinations of those given in ref. ||T5[ (again 
we have discarded subleading terms which are the product of two traces). Note that 
the last term is not immediately obtained by inserting A_ within the trace in ci'^^; in 
order to obtain ( pl|) , one has to use Df^T. = -S (-D^St) S and 1 = S St ^ S A_ St. 
In the case of C^^\Qii) it was sufficient to calculate the K ^ tt transition induced by 
Qii to determine the unique coefficient. Here we have to find three coefficients instead. 

The i^T — > vr transitions are given by: 



2i 



)e,h = -p + "^d) k G 



(4) 



(4.13) 
(4.14) 



which show that (Qe) and G^^^ (Qe) cannot be distinguished at this level (to consider 
other off-shell meson-to- meson transitions does not solve the problem). Thus, from the 
calculation of X — > vr transitions (see Fig. 3 for the relevant diagrams) we can at 
most determine G^^\Qq) and the sum 

[G^EiQe) + Gn^Qe)]- We therefore have also 
to consider the — > 27r amplitudes at the quark level, and match them with the chiral 
lagrangian results: 

^(m,-m,)m2Gg\ (4.15) 
^ {nis + nid) (m|- - ml) 

md)ml 



vr ' TT 



A{K^ 



+ {rriu 
1^ 



G 



(4) 
H ' 



{rud - niu) ml + {m\ - ml) 



G 



(4) 
K 



(4.16) 
(4.17) 



From eq. ( 4.16| ) it appears that in order to determine G^^ {Qq) by calculating quark 
loops for K 2tt, we may keep the m^m^ terms only. Moreover, if the coefficients 



G 



(4) 



E H k(Q^) same order of magnitude, the term will be the most impor- 

tant one. 

(4) 



To determine the coefficients GK 



EHxiQ^)^ we calculate the K ^ it and K ^ 2tt 
amplitudes due to Qq within the xQM. For the i^T — > vr transitions, in the leading 
factorizable limit, we obtain: 

(7r+|g6|i^+) 



2(vr+|ll75d|0)(0|s75n|K+) 



{0\dd\0) + (0|ss|0)l {7T+\sd\K+) 



(4.18) 
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and for X — > 27r: 

(7r+^-|Q6|i^°) 



2 (7r+|¥75d|0)(7r-|su|/^°) - 2(^+7r-|(id|0)(0|s75(i|-fC° 



+2 



(0|ss|0) - {0\dd\0) (7r+7r~|s75d|i^°). (4.19) 

These equatfons contain some building blocks which we calculate within the xQM. 
The quark condensate (Fig. 2) with mass insertions included is given by 



im) = {qqf^ + ^ims+mu)[{qqf^ + M//W] + {qq)G . 



(4.20) 



where the last term represents the two-gluon condensate contribution (see eq. ( |3.10D ). 
As for the diagrams in Fig. 4, in the naive dimensional regularization (NDR) scheme 
(anti-commuting 75 in d 7^ 4), we obtain: 



(0|s75u|i^+(fe))NDR 



iV2 



k 



(0) 



2 ./TT 



2M 



{7r+ip+)\sd\K+{k)) 



NDR 



3M 2 , ^/ 
/2 ^ 2A2 2M 



(4.21) 



/4 



X 

- I /+ - 6 ^ 

.(0) 



M2 



q^ + q- P 



(4.22) 



where q = k — and P = A; + while = /j^ // = 1 can be identified with the 
vector form factor at zero momentum transfer q. Similarly, in the 't Hooft-Veltman 
(HV) scheme (commuting 75 in d 7^ 4), we find: 



{0\sj,u\K+{k))^, 



scriptsize HV 



(0|s757x|i^+(A;))NDR 



iV2f 



+12 {nis + ruu 



M3 / 



6M2 



{7T-^{p+)\sd\K+ik)) 



HV 



{7r+{p+)\sd\K+ik)) 



NDR 



24 



M3 
A2" 

X 



In eqs. ( [4.2lD -( ^.24| ) we have introduced the chiral symmetry-breaking scale 
~6 



A. 



27r 




7 = 0.83-1.0 GeV, 



(4.23) 
(4.24) 

(4.25) 
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where the given range corresponds to assigning to / the numerical values of /t^ and fx, 
as indication of S'C/(3)-breaking effects. Notice that identifying = / in eq. (2.5) 
leads to / ~ ^/N^, which makes independent of N^. As a consequence, the 
dependence of eqs. ( [4.21| )-( p^ ) resides entirely in / and (qq) ~ N^.. 

The terms containing current quark masses are obtained from mass-insertions due 
to Cmass in eq. (0) in the various quark loop diagrams in Fig. 4. We have discarded 
terms proportional to rrici because they are not needed in determining G^ehk^Q^)- 
The matrix element (7r^(p+)|ll75d|0) is obtained from {{)\'s^^u\K^ [k)) by obvious sub- 
stitutions, while (7r^(p+)7r^(p_)|dd|0) (where p+ + p_ = k) can be derived from 
eq. (|]22|) and ( g^ by crossing {k —P-)- Finally, (vr+vr {sj^dlK ) is zero to 
this order p7| , Notice that these expressions are more general than those given 
by other authors |2^, ^] , because our results are not based on the divergence of the 
on-shell current. 

By using these matrix elements, we find that the constant term for the Qg— induced 
-fC — > TT transition vanishes as it should according to chiral symmetry p^, 29, |30|]. This 



cancellation is exact in NDR; in the HV, it holds only up to terms of order 1/A^, a 
warning about the relevance of higher-order loop effects. A typical example of such 
a high order contribution would be any diagram in Fig. 3 with an extra meson line 
connecting the two quark loops. For consistency, we will always drop terms of order 
1/A^ or higher in all numerical estimates. 

In the NDR scheme, neglecting current quark masses which would amount to a 
correction < 10%, the LO coefficient for the Qq matrix element is therefore 



Gf (%) 




while the result in the HV scheme is 
2/2 



Gg'^ {Qe 



M 



im) 



(0) 



/4 



M2 
X 



(gg)(o)A2 



1 + 4 



Then the relation 

.(2) 



G's'iQe) 



(4.26) 



(4.27) 



(4.28) 



defines the NLO chiral parameter L5 in agreement with refs. 0, In passing, let 
us remark that the leading non-zero matrix element of Qq is next-to-leading in the 
framework of the chiral perturbation theory since it is proportional to the coefficient 
L5 of Tr [(St Mg + Ml S) D^^StD^S] in ^i^^ . 
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In the xQM, to order M^/A^, we find 
^r°^ = -^llLf/+-6^'l , (4.29) 



X 



SM{qq) A2 
as in ref. [Bl, and 



L^^- f 2^\l2^^\ (4 30) 

botli sensitive functions of M. Let us remark that for 

eq. ( [4.28 ) used in eq. ( ^ ) gives the usual l/A'c matrix element Q. In that scheme we 
therefore have 

L5 = 1.4x10"^ for A^ = 0.83GeV 

L5 = 1.0x10"^ for A^ = 1.0GeV. (4.32) 

A scale-dependent value for L5 can be extracted in xPT from the ratio fx/fn where 
one finds |31] 

L^inip) = (1.4 ± 0.5) X 10-3 , (4.33) 

which we use for comparison with eqs. (4.29)-( 4.3lD . 

In Table 2 we have collected the results of the xQM predictions for L5 in the two 
schemes as we vary M. A constraint of approximately 3a from the central value in 
eq. ( |4.33| ) gives 

140 MeV ^ M ^ 340 MeV , (4.34) 

where we have taken for the scale dependent quark condensate the expression in 
eq. (5.12), and assumed a perturbative running for /i > 1 GeV. Recall however that for 
M > 300 MeV M^/A^ corrections become relevant and may affect the smaller values 
of L5 listed in Table 2. 

Moving on to the NLO contributions, by using eqs. (4.15)-(4.17) we find in the 
NDR scheme: 

g2'«3.) = 3Cafe>'«)^ + ^^), (4.36) 
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whereas in the HV scheme: 



/2 _ 



M{qq)(^) 



1 + 12' 



A2 

X 



(4.37) 



1 + 12 



(g-<7)(0)A4") 




6M(gg)(0) 



+ 12' 



A2 

X 



(4.38) 



The most hnportant contribution to K ^ 2tt comes from the third coefficient 
G^^ {Qo), since the corresponding amphtude is enhanced by a factor m\/m\ compared 
to the other two (see eqs. ( 4.15| )~( 4.17] )). In the NDR scheme we find: 



G^H (Qg 



77n\(0) 



(4.39) 



whereas the HV result is: 



Gff^ {Qe 



-7Ce{qq) 



(0) 



f 156 M^f 



7 



7M{qq)(^') 



(gg)(o)A4 
/+- 




1 + 12 



A2' 

X 



(4.40) 



Recently an NLO analysis of the hadronic matrix elements Qi-io has appeared |3^, 
which makes use of effective scalar-meson exchange. The result of ref. for the 
operator Qq does not contain terms corresponding to (with two derivatives and 
current masses squared), which contribute to eqs. (4.15)-( fl.l7D . 

We see that G^^^jj j^{Qq) are formally suppressed by terms proportional to either 
f/Al or f/{M{qq)^^^) with respect to G(^)(Qii). However, because of the numerical 
factors in front, the NLO G^'^\Qq) coefficients turn out to be numerically of the same 
size as G^^^Qu). Since the amplitude due to £-^^\Qn) is effectively suppressed by 
a factor m^/m?^ with respect to the amplitudes obtained from C'^'^\Qq) and Ch, the 
contribution of Qn to e' /e turns out to be generally more than one order of magnitude 
smaller than these NLO corrections (see Table 3). The same can be said of other 
subleading operators such as Qa, and Qj (see eq. ( jST^ ) below for the definitions, 
and Tables 7-9). 

In Table 3 we have reported the weights of the NLO matrix elements of and Qii 
relative to the LO amplitude, as computed in the xQM. These estimates are robust 
as they do not sensibly depend on the detailed values of Aqcd and A^. 



15 



5 Scheme dependence of (Q6)8 and {Q^)27 



The Wilson coefficients computed by means of the renormalization group equa- 
tions depend, at the NLO, on whether the NDR or the HV prescriptions are used in 
treating 75 in d dimension. At the same time, the matrix elements of the relevant 
operators do not have any scheme dependence, at least if computed by 'i/Nc tech- 
niques. As a consequence, there remains a 75-scheme dependence which increases the 
uncertainties in the final result. The discrepancy between the HV and NDR result 
goes from 30% {Kqcd = 200 MeV) to 40% (Aqcd = 400 MeV) at mt = 130 MeV, 
becoming much worse in the "superweak" regime at mj = 170 GeV, where it goes from 
70% (Aqcd = 200 MeV) to 80% {Kqcd = 400 MeV) §. 

A potentially nice feature of the xQM is that it makes it possible to compute also the 
matrix elements in both schemes. In eqs. ( [4.26| )-( ^.27| ) we have computed the matrix 
element for the operators Qq to the leading order in the chiral expansion. The next 
most relevant contribution comes from the electroweak penguin operator Qg (eq. ( |1.3D ). 
This operator can be written as 



where 



-2ed 



2/3, so as to obtain 



(7r+7r-|Q8|i^°) 



1 



{7T+Tr-\QQ\K^ 



and, after a Fierz rearrangement and factorization, 

(7r+^-|Q8|-f^°)27 = 3 (7r+|n75tf|0)(7r-|sn|K°) , 
(7rV|Q8|i^°)27 = 0. 
The leading order amplitude corresponding to eq. ( |5.3| ) is therefore given by 



^(i^°^vr+7r-; Q8)27 
in the NDR scheme, and 

in the HV scheme. In eqs. (| 



3\/2 



3^2 



(0) 



(5.1) 
(5.2) 

(5.3) 
(5.4) 

(5.5) 



1 + 12 



(5.6) 



|) we have kept only the leading momentum- 



independent terms. The momentum-dependent corrections are at the 10% level. 
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In a toy model for e' je in which we use the full mixing to determine the Wilson 
coefficients, but keep only the contribution of the two leading operators Qe (for the 
1 = amplitude) and Qg (for / = 2), we can compare the scheme dependence in the 
l/A'^c approach — where it only appears in the Wilson coefficients — and in the xQM, for 
different choices of M. 

We have taken the Wilson coefficients at ;U = 1 GeV (/x = 1.4 GeV) for rrit = 
170 GeV and multiplied them by the corresponding matrix elements, according to 
eq. ( |6.5D in the next section. We have used 



I (^+7r-|Q,|i^0) 



+ 



I (vtVIQ.IETO) 



{7r+7T-\Qi\K' 



{7r+7r~\QQ\K' 



^{7T+7r~\Qs\K^)27. 



{7t\^\Q,\K^) 



(27r, / = 0|Qi|K°) = 

{2n, I = 2\Qi\K'>) = 
In our case 

(2^, / = 0|Q6|^°) = 

{2tt, I = 2\Q8\K'^) = 
Consistently to our approximation, we neglect the / = contribution of Q 

{2tt, / = 0| Qs |i^°> = - I (2vr, / = 0| Qe \K° 



+ 



^(vr+vr-IQsli^O; 



27 



(5.7) 
(5.8) 

(5.9) 
(5.10) 

,(5.11) 



since it is suppressed in e'/e by a factor 1/20 (A/ = 1/2 rule) with respect to 
the 1 = 2 contribution of eq. ( [5.10| ). Notice that, at the zeroth order in momentum 
expansion, (vr+vr^lQel-^'^) = and, as a consequence, the last term in eq. ( ^.11 ) is 
leading. 

The matrix elements in the xQM are computed using the scale-dependent quark 
condensate 



(99) (/^) 



(5.12) 



Following ref. |], we take m^(1.4 GeV) = 150 MeV and mrf(1.4 GeV) = 8.0 MeV. The 
values at other scales are obtained using the NLO evolution of QCD. In particular, for 
Aqcd = 300 MeV one has m^(1.0 GeV) = 176 MeV and mrf(1.0 GeV) = 9.4 MeV. 

Tables 4 and 5 summarize our findings for different values of Aq^^. In using 
eq. (4.27) and eq. ( p. 61 ) we have dropped all terms of order 1/A^ for consistency and 
limited ourselves to values of M < 220 MeV. As explained in the beginning, we cannot 
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trust our results beyond order 1/A^ because we have neglected higher- loop corrections 
with meson exchange that give 1/A^ contributions. 

In the range considered, the 75-scheme dependence is indeed dramatically reduced. 
For instance, taking Aq^^ = 300 MeV, A — defined as the difference between the 
HV and the NDR results divided by the HV one— is below 10% for M ~ 120 - 150 
MeV (// = 1.4 GeV) and for M ~ 120 - 170 MeV (/i = 1 GeV), when A^^ = 830 MeV; 
similarly, A is below 10% for M ~ 120-180 MeV (/x = 1.4 GeV) and for M ~ 150-190 
MeV {fi = 1 GeV), when A^ = 1 GeV. 

Another, and more restrictive, reading of the same tables puts together all data 
for different Aqqd^s. In this case, stability is achieved for M ~ 120 — 150 MeV (/i = 
1.4 GeV) and for M ~ 150 MeV {fi = 1 GeV), when A^ = 830 MeV; similarly, A is 
below 10% for M ~ 120 - 170 MeV {fi = 1.4 GeV) and for M ~ 190 MeV = 1 GeV), 
when A;^ = 1 GeV. 

To reach values of M larger than its central value and closer to it, we would have 
to include higher-loop corrections. The preliminary nature of our analysis needs hardly 
be stressed as only two out of eleven operators have been considered. Even though Qq 
and Qs induce the most relevant contributions, we expect that a complete estimate of 
e'/e would result in a more stable range of values of M for which the 75-scheme and /j, 
dependences are reduced. This will also give us confidence on the size of e'/e that is 
obtained. 



6 "NLO" STUDY OF e'/e with dipole operators 

In this section we discuss the impact of the dipole gluon penguin Qn on present 
estimates of e'/e. The contribution of Qii is given by setting {qq)G at the rather 
conservative value of — (275MeV)^ (see Table 1). 

Since a satisfactory calculation of e'/e in the xQM is missing beyond the leading 
factorizable order (the study is under way), we resort to the 1/Nc analysis of ref. 
for the ten traditional operators. We present our results in tables that show a detailed 
anatomy of the contributions of the different operators, and can be used as reference 
for future developments. 

Because of the many ingredients involved in the calculation of e'/e, it is useful 
to briefly recall the theoretical inputs used. The effective lagrangian for |A5| = 1 
transitions can be written, for /x < rric, as ||2| 

^AS=i = ^A„^U(/x) +ryi(/i) Qi(/x) = ^Ci(/i)Qi(/x) . (6.1) 

V2 ,• L J 
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In the previous equation, Aj = VidV*,, where V is the Kobayashi-Maskawa (KM) matrix, 
and T = —Xt/Xu- The Wilson coefficients 21^2 (^) run from rriw to rric via the corre- 
sponding 2x2 sub-block of the 10 x 10 anomalous-dimension matrices, while Zi{fi) = 
for i = 3 — 10. From fi = rric down, as the charm-induced penguins come into play, all 
Zi{ii) evolve, given the proper matching conditions, with the full anomalous-dimension 
matrices. The Wilson coefficients {Viifj) = Vi{fj,) — Zi{fi)) arise at m]y due to 

integration of the W and top quark fields. They coincide with Zi(fi) for i = 1,2, the 
information about the top quark being encoded in the z = 3 — 10 components. 

The list of the effective operators Qi {i = 1 — 10) is reported in refs. |§, |3|, whose 
notation we follow closely and where the reader may find a complete discussion of the 
basic tools. For convenience we report here the ten operators usually considered 

Ql = (Sa^i/3)v-A (%^a)v-A ' 

Q2 = (s^i)v-A M)v-A ' 

(53,5 = (sd)v-A Eg {m)y^A > 

(6.2) 

Q4,6 = {Sadf3)y^j^J2q{Q^Qa)vTA, 

Q7,9 = i (srf)v-A Eg eg (Qg)v±A > 

Qs.lO = i (SaC^/3)v-A Eg eg(^/39a)v±A , 

where a, (3 denote color indices (a,/3 = l,...,A'c) and eg are quark charges. Color 
indices for the color singlet operators are omitted. {V it A) refer to 7^(1 ± 75). We 
recall that Qi^2 stand for the 1^-induced current-current operators, Q3-6 for the QCD 
penguin operators and Qt-io for the electroweak penguin (and box) ones. 

Not all the operators in eq. ( |6.2| ) are independent. For /i < mc, having integrated 
out the charm quark, we have 

Qi = Q3 + Q2 — Ql , 

Qg = (3Qi-Q3)/2, (6.3) 
Qio = Q2 + {Ql - Q3) /2 . 

Note that these relations hold in the HV scheme, but they may receive additional 
contributions in other schemes since Fierz transformations have been used in obtaining 
them. 

Together with this basis, which closes under QCD and QED renormalization, one 
should a-priori consider two additional dimension-five operators: Qn (eq. (|1.1| )) and 

Qi2^^m,sa^,F'^'^(l-75)d, (6.4) 
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which account for the magnetic and electric dipole part of, respectively, the QCD and 
electromagnetic penguin operators. In eq. (6^) = —1/3 is the charge of the down 
quarks. 

Actually, in Q we argued that the hadronic matrix element of the electromagnetic 
operator ( |6.4D is negligible and, accordingly, even if we keep the operator in our basis 
for the Wilson coefficients, we will put its contribution to zero in the end. 

Since lm{Xu) = according to the standard conventions, the short-distance com- 
ponent of e' /e is determined by the Wilson coefficients yi. Following the approach of 
ref. Q, yi(/i) = 1/2 (/^) = and the effect of (5i,2 appears only through the linearly- 
dependent operators Q4 9 lo- 



e 
e 

where 



The lagrangian in eq. (6.1) yields 
Im At 



10 



-4 



1.7 X 10-4 



p(l/2) _ p(3/2) 



p{l/2) 
p(3/2) 



r ^yi(27r, / = 0| 
-5]yi(27r, I = 2\Q,\K'^) 



[I -a 



(6.5) 

(6.6) 
(6.7) 



We take, as input values for the relevant quantities, the central values given in ap- 
pendix C of ref. IQ]. This allows us to reproduce, in the ten-operator case, the central 
values of the results given in appendix B of ref. |^]. In particular, we take 
ujGf 



1.7- 



594 GeV" 



1/22.2 



i^T^-i-T^/ — 0.25 ] 



(6.8) 



2|e|Re^o 

Im At is determined from the experimental value of e as an interpolating function of 
rrit- Its central value, given the KM phase 6km in the first or second quadrant, is given 
by 



and 



ImAt ~ 2.77 X 10 



ImAt ~ 2.19 X 10" 



-4 



-0.365 



-0.47 



(first quadrant) , 



(second quadrant) 



(6.9) 



(6.10) 



where xt 

The value of the Wilson coefficients yn and yi2 at the hadronic scale of 1 GeV can 
be found by means of the renormalization group equations. Denoting generically the 
vector of Wilson coefficients by C(/^), its scale dependence is governed by 



dg. 



C 



m 



W 2 

g ,a 



l^{g\a)C 



m 



w 2 
■^,g ,a 



(6.11) 
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where P{g) is the QCD beta function and a the electromagnetic coupling (the running 
of a is being neglected). At the NLO we have 

7b^") = S(^^°' + £tf')+^(7f + ^7£') , (6.12) 

where 7!°^ and 76^^^ govern the leading QCD and the electromagnetic running respec- 
tively. The anomalous-dimension matrices labelled with (1) refer to the NLO running 
(O(a^) effects are neglected). 

In order to include all available NLO effects in the evaluation of e' /e, we follow the 
analysis described in ref. The NLO 10 x 10 mixing matrices for the operators Qi-io 
can be found in refs. ||2|, ^. Concerning the dipole operators, the leading-order matrix 
of the strong anomalous dimensions of Qii and Q12 and their QCD-induced mixing 
with Qi_6 can be borrowed from the existing calculations for the b ^ decay |33] 
(recent discussions are given in ref. |]l2[). In fact, by replacing s ^ b and d — > s in 
eqs. (|6.2| )- (|6.4D we obtain the operator basis, which should be considered for a complete 
NLO analysis of b ^ 57. 

While the 10 x 10 part of the anomalous dimension matrices (6.12) is identical to 
that used in refs. |2|, |^, two extra columns and rows have to be added to represent the 
mixing of the first ten operators with the two new ones, which takes place first at the 
two-loop level. We have taken for all two-loop anomalous dimensions the HV scheme 
results @, In this way, no finite additional contributions to the renormalization of 
yii and yi2 arise at the various quark thresholds (for a discussion see Misiak in ref. ||3^ 
and ref. 0]). The explicit expression of 7J is reported for instance in ref. We 
just recall that whereas the evolution of the dipole Wilson coefficients Cn and C12 is 
substantially affected by the mixings with Qi-io, the Wilson coefficients Ci_io remain 
unaffected by the presence of the dipole operators (Qi-io close under QCD and QED 
renormalization). The two-loop mixings of Qn and Q12 with the electroweak penguins 
(Qj-io), not yet given in the literature, can be easily derived from those with the gluon 
penguins (Qs-g)- We have verified that their effect on the running of Cn and C12 is 
negligible (< 1%) and therefore they can be safely set to zero (and, by extension, they 
can also be set equal to zero in 6 — > S7). 

The complete NLO analysis would require computing, among other things, the 
three- loop mixings between the dipole operators and the first ten (quite a task!). The 
lack of knowledge on these entries introduces an uncertainty on the dipole Wilson 
coefficients which can be as large as 50% (see the analogous discussion for the 6 — > 57 
inclusive decay in ref. |jl2|). However, for what concerns e' /e this uncertainty is diluted 
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over many contributions, and it is certainly not as relevant as our ignorance of the 
hadronic matrix elements. The results presented here are obtained by adding two rows 
and two columns of zeros to the 10 x 10 electromagnetic and NLO anomalous dimension 
matrices, which can be found in refs. As a consequence, the contribution of Qn 

to e'/e takes into account only the "LO" (two-loop) QCD effects. 

The expressions for the initial Wilson coefficients vi_io iTnyv) can be found for 
instance in ref. [Q]. For what concerns the new coefficients f 11,12 (nT-w) we have 

vn{mw) = -E'{mf/mly) and vuimw) = -D'{m^/m^r)/ed , (6.13) 

where 

E'(x) = — -Inx (6.14 

^ ^ 2(1 -x)4 4(1 -x)3 ^ ' 

x^(2 - 3x) , 8x3 ^ 5^2 _ 7^ 

D (x) = ; 71-lnX ; 7^ — . (6.15) 

^ ^ 2(1 -x)^ 12(1 -x)3 ^ ' 

In Table 6 we report the HV results we have obtained for 21-12 (a* = 1 GeV) and 
2/3-12 (/^ = 1 GeV) (recall that yi,2 (/f^) = 0) compared with their initial values, for 

= 170 GeV and for A^J.^ = 200, 300 and 400 MeV. We fully agree on the values of 
the renormalized coefficients for the first ten operators with ref. [||. 

Let us remark that the effect of operator mixing induces a renormalization on ^11,12 
(= yii,i2 + -211,12) which is a factor of 4-5 larger than that induced by multiplicative 
running alone (which roughly reduces by a factor of two the initial Wilson coefficients) . 

In order to discuss the effect of the various operators in determining the size of 
e'/e, we need a consistent estimate of the relevant hadronic matrix elements. For the 
operators Qi_io, we follow the strategy of ref. |^] where the various matrix elements 
are evaluated by means of the 1/Nc expansion and soft-meson methods. Overall co- 
efficients and B^^^"^^ {i = 1 — 10) parametrize our level of ignorance of their 
normalization scale, scale dependence and the approximation inherent in the method. 
The matrix elements of Qi and Q2 can however be determined phenomenologically 
from the experimental values of ReAg and JieA2, so as to reproduce the AI = 1/2 
rule. In particular, in ref. ^ it is found that at /i = nic, B^j^^Q « 6.3 in the HV 
scheme, which is about three times larger that the 1/Nc result. Related to this coeffi- 
cient is the value of B^j^j^Q, which we find equal to 20.2, 13.5 and 8.1 for a[J^^ = 200, 
300 and 400 GeV, respectively. Correspondingly, ^(^/r^^o = 0-^5, 0.46 and 0.48. The 
large deviations from unity of these effective coefficients gives us a gauge of our lack 
of understanding of the AI = 1/2 rule within the 1/Nc expansion (better results are 
achieved within the xQM P], leaving deviations at most of a factor 2-3). 
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Further relations among other coefficients are advocated in ref. depending on 
the relevance and the role of the various operators, so as to reduce, in the ten-operator 
case, the description of the hadronic sector to two effective parameters: B^^"^^ and 
B^^'^\ whose leading 1/Nc value is 1. 

The inclusion of Qu and Q12 requires three additional effective parameters: By/'^\ 
B^j'^'^ and Bfj'^\ For Qn we use our result (|3.1lD , where for {qq)G we take a typical 
value from Table 1, namely —(275 MeV)^. 

Since the determination of Bi and B2 is best achieved at // = rric all the 
hadronic matrix elements are assumed to be evaluated at that scale and renormalized 
down to 1 GeV via their anomalous-dimension matrix. We proceed analogously, by 
setting B^i'^^ = 1 and, as we neglect Q12, b'^J'^'^ = ByJ"^^ = at = rric and using the 
12 X 12 QCD and electromagnetic evolution matrices to evolve all the hadronic matrix 
elements to the 1 GeV scale. 

Since the anomalous-dimension matrices, which govern the evolution of the hadronic 
matrix elements, are the transpose of those evolving the Wilson coefficients, we now 
find that the presence of Qii,i2 affects, from /i = rric down, the renormalization of the 
first ten operators. On the other hand, the evolution of (5ii,i2 is determined solely by 
their 2x2 anomalous-dimension matrix, which implies that the matrix element of Q12 
remains vanishing. 

As a consequence of the previous remarks, our results for the individual contribu- 
tions of the operators Qi-io to e' /e may differ slightly from those reported in ref. [^. 
We have however checked that, in the ten-operator case, we reproduce their NLO results 
exactly. 

Tables 7, 8 and 9 show the contributions to e' /e of each operator, for different 
choices of J^q(J£, and nit, in the HV scheme. The first ten contributions are also 
partially grouped in a "positive" gluonic component versus a "negative" electroweak 
component, which shows the "superweak" behavior of e' /e within the standard model 
as the top mass increases. The total effect in the extended operator basis is given in 
the last row. 

We find these tables a useful way of displaying the currently available theoretical 
information on e'/e. In particular we observe that for mt ^ 170 GeV the size of e' /e 
becomes comparable to the Qn contribution alone, signalling the relevance of NLO 
contributions to the hadronic matrix elements. 
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A Feynman rules 



In this appendix we collect some formulas that are useful in computing the hadronic 
matrix elements. 

The free propagator for the constituent quark is given by 



Soip) 



(A.l) 



where ^ = 'j ■ p- The same propagator in the external gluon field (fixed-point gauge) 
is 11: 

Slip) = - ^ G%t^ ..2^2 ' (A-2) 

where 



(A.3) 



and cj^i, = {i/2)[y^, j^]. In order to compute the gluon condensate, the quark propa- 
gator in two external gluon fields is needed: 

where 



(p2 _ JV^2)5 

= 7a(i' + Af)7/3(^ + Af)7^(i^ + M)7,(^ + M) . 



Other useful formulas are: 
6 ^ vr 



and 



cj'^V^, = 12 I : 



(A.5) 



(A.6) 



(A.7) 



The relevant meson-quark interactions are derived from the lagrangian in eq. (|27 
which we write here as 

M „ „ M 



^xQM = -M^(7 + 2iyg75n(? + 2^gn2g + 0(l//3) 



where 



2^ V2 





7r+ 














^8 



(A. 



(A.9) 
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and 

which yields 

^75 Hq = -j= (uj^uTT^ — d^^diiP + d^^sK^ + u'y^d'K'^ + • • • j . (-^-H) 

The relevant Feynman rules are therefore: 

n - M^/2 
d75S - coupling: — 75 (A. 12) 

M\/2 

wy^s — coupling: j — 75 (A. 13) 

7r° dj^d - coupling: '^~J''^^ i-^-^^) 

TT^ uj^u — coupling: — ^75 (A. 15) 

+ - , ,. MV2 

'"75" ~ coupling: — 75 (A. 16) 

7r° ds - coupling: "^/^ ^'^''^^^ 

7r+ us - coupling: (A.18) 

TT^ US — coupling: +i (A. 19) 

/ v2 

, - M 
TT ds- coupling: . (A.20) 

All meson fields are entering the vertex. The same rules hold for the conjugate cou- 
plings. 
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IVl I iV±C V J 


1 9(1 


140 




1 sn 

-LOU 




990 


{qqff (MeV) 


-315 


-299 


-286 


-275 


-265 


-257 


M (MeV) 


240 


260 


280 


300 


320 


340 


{qq)T (MeV) 


-250 


-243 


-237 


-232 


-227 


-222 



Table 1: Values of {qq)G = ~J2M i^^^) ^ function of M. The central value 
{^GG) = (460 MeV)'^ has been used. The scaling of the given entries for other values 
of the gluon condensate is straightforward. 





is X 10^ 










L5 X 10^ 






(/X = 1.4GcV, {qqy/^ = 


-273 McV) 




in = 1.0 GeV, {qqy/^ = 


-259 MeV) 


Ax (GeV) 


0.83 


1.0 




Ax (GeV) 


0.83 


1.0 


M (MeV) 


NDR 


HV 


NDR 


HV 




M (MeV) 


NDR 


HV 


NDR 


HV 


120 


3.4 


3.7 


3.6 


3.7 




120 


4.0 


4.3 


4.2 


4.4 


140 


2.8 


3.1 


2.9 


3.2 




140 


3.2 


3.6 


3.5 


3.7 


160 


2.3 


2.6 


2.5 


2.7 




160 


2.7 


3.0 


2.9 


3.2 


180 


1.9 


2.2 


2.1 


2.4 




180 


2.2 


2.6 


2.5 


2.7 


200 


1.5 


1.9 


1.8 


2.1 




200 


1.8 


2.2 


2.1 


2.4 


220 


1.2 


1.6 


1.5 


1.8 




220 


1.4 


1.9 


1.8 


2.1 


240 


1.0 


1.4 


1.3 


1.6 




240 


1.1 


1.6 


1.5 


1.8 


260 


0.7 


1.2 


1.1 


1.4 




260 


0.9 


1.4 


1.3 


1.6 


280 


0.5 


1.0 


0.9 


1.2 




280 


0.6 


1.1 


1.0 


1.4 


300 


0.3 


0.8 


0.7 


1.1 




300 


0.4 


0.9 


0.8 


1.2 


320 


0.2 


0.7 


0.6 


0.9 




320 


0.2 


0.7 


0.7 


1.0 


340 


0.0 


0.5 


0.4 


0.8 




340 


0.0 


0.5 


0.5 


0.9 



Table 2: Table of L^^j,) obtained in the xQM for different values of M. The results 
are shown for A^l^. = 300 MeV. 
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M (MeV) 


120 


140 


160 


180 


200 


220 




{Q6)ndR 




1 % 


1 % 


1 70 


1 % 


2 % 


2 % 








1 % 


1 % 


1 % 


1 % 


2 % 


2 % 




/n \NLO 
W6/NDR 

{Q6)nDR 




14 % 


15 % 


17 % 


20 % 


22 % 


25 % 








14 % 


16 % 


18 % 


20 % 


22 % 


25 % 



Table 3: Typical weights of NLO contributions to the K tt'^tt amplitude in the 
XQM. The results arc shown for ^ = 1 GeV, a[^^^ = 300 MeV, = 1.0 GeV and 
mt = 170 GeV, and include the NLO Wilson coefficients ("LO" for Qn). 
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"e'/e " xlO^ (m = 1-4 GeV, = 0.83 GeV) 




200 McV 


300 MeV 


400 McV 




NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 




8.7 


7.3 


-19 % 


11 


9.0 


-22 % 


14 


11 


-26 % 


M(MeV) 


NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 


120 


21 


21 


-1 % 


27 


26 


-6% 


34 


32 


-6 % 


140 


15 


16 


4 % 


19 


20 


2 % 


24 


24 


-0 % 


160 


10 


12 


14 % 


13 


15 


12 % 


16 


18 


9 % 


180 


6.5 


9.2 


30 % 


8.1 


11 


28 % 


10 


14 


26 % 


200 


3.2 


7.0 


54 % 


4.0 


8.6 


54% 


5.0 


10 


52 % 


220 


0.4 


5.4 


92 % 


0.4 


6.5 


94% 


0.5 


7.9 


94% 





"e 


/e"xlO'* (/X = 1.0 GeV, Ax 


= 0.83 GeV) 








200 MeV 


300 MeV 


400 MeV 




NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 




8.0 


6.5 


-24 % 


10 


7.6 




-33 % 


13 


8.7 


-46 % 


M(MeV) 


NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 


120 


23 


22 


-4 % 


31 


28 




-10 % 


43 


37 


-18 % 


140 


17 


17 


1 % 


23 


22 




-5 % 


32 


28 


-14 % 


160 


12 


13 


11 % 


16 


17 




4% 


23 


22 


-6 % 


180 


7.5 


10 


25 % 


10 


13 




17% 


16 


17 


7 % 


200 


4.1 


7.8 


47 % 


6.3 


10 




38 % 


11 


14 


26 % 


220 


1.1 


6.1 


81 % 


2.4 


8.1 




70 % 


5.1 


11 


54% 



Table 4: Toy model for the 75-sclieme dependence of e' je. The results are shown for 
mt = 170 GeV, A-^ = 0.83 GeV, and 5km in the first quadrant. The two tables refer 
to the choice ^ = rric =1.4 GeV and 1.0 GeV respectively for the renormalization 
scale. The xQM results are compared with the 1/Nc predictions evaluated according 
to ref . [2] . The parameter A is defined as the (HV — NDR) /HV combination of the 
entries. 
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V/e " X 10^ (m = 1-4 GeV, = 1.0 GeV) 




200 MeV 


300 MeV 


400 MeV 




NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 




4.4 


3.3 


-33 % 


5.5 


4.0 


-37% 


7.0 


4.9 


-43 % 


M(MeV) 


NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 


120 


23 


21 


-6 % 


29 


26 


-8% 


36 


32 


-11 % 


140 


17 


16 


-3 % 


21 


20 


-5% 


27 


25 


-8 % 


160 


12 


13 


2 % 


15 


15 


-0 % 


19 


19 


-3 % 


180 


8.7 


9.7 


11 % 


11 


12 


9 % 


14 


14 


5 % 


200 


5.6 


7.5 


24 % 


7.1 


9.1 


23 % 


8.9 


11 


20 % 


220 


3.1 


5.7 


46 % 


3.8 


7.0 


46 % 


4.8 


8.4 


43 % 





V/e " X lO'' 


[jjL = 1.0 GcV, Ax = 1.0 GeV) 






^(4) 


200 MeV 


300 MeV 


400 MeV 




NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 




4.0 


2.7 


-45 % 


5.0 


3.1 


-60 % 


6.4 


3.5 


-84 % 


M(MeV) 


NDR 


HV 


A 


NDR 


HV 


A 


NDR 


HV 


A 


120 


25 


23 


-9 % 


33 


29 


-15% 


46 


38 


-23 % 


140 


18 


17 


-6 % 


25 


22 


-12 % 


35 


29 


-21 % 


160 


14 


13 


-1 % 


19 


17 


-8 % 


27 


23 


-17% 


180 


9.8 


10 


7 % 


14 


14 


-1% 


20 


18 


-11 % 


200 


6.6 


8.2 


19 % 


9.7 


11 


11 % 


15 


15 


-1 % 


220 


3.9 


6.5 


39 % 


6.1 


8.6 


29 % 


10 


12 


15% 



Table 5: Same as Table 2, with = 1 GeV. 
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ZfUU ivit; V 


OUU iVlc V 


4nn MpV 


Zl 


(U.Uolj —0.477 


/A AOO\ A £?A/? 

(0.033) — 0.6U6 


(0.035) —0.774 


Z2 


(0.98oj 1.250 


/A AOO\ 1 Ad 

(0.988) 1.346 


(0.987) 1.472 


Z3 


U.UUD 


U.UUo 


n ni 1^ 

U.UiO 


Z4 


— n ni 1 


— n ni R 


— n n'^n 

u.uou 


Z5 


u.uuo 


n nnzi 


n nnfi 

u.uuu 


^6 


— n ni n 


— n ni fi 

u.uiu 


— n n9fi 




— n DfTi 

\J .\J\JO 


— n nn3 


— n nn9 


z^/ a 




01 1 


01 ^ 




—0 001 


004 


010 


Zio/ Oi 


— n nn? 

U.UU / 


— fl ni 1 

U.U J- 1 


— n ni 7 

U.U-L / 


Zll 


U.UOtJ 


— n (\AA 


— n n^fi 

u.uou 


Z\2 


u.o^^ 




n fiQn 

U.UiyU 


2/3 


(O.OUlj U.UzD 


/A Am \ A AO/( 

(O.OUlj 0.034 


(U.UOUj U.044 




(0.001) -0.01c 


(0.001) -0.05(3 


(0.001) -0.066 




(0.000) 0.013 


(0.000) 0.015 


(0.000) 0.018 


ye 


(0.001) -0.066 


(0.001) -0.088 


(0.001) -0.120 


Vila 


(0.151) -0.031 


(0.151) -0.029 


(0.151) -0.027 




(0.000) 0.126 


(0.000) 0.172 


(0.000) 0.240 


yg/a 


(-1.094) -1.541 


(-1.094) -1.632 


(-1.094) -1.759 


yio/a 


(0.000) 0.560 


(0.000) 0.703 


(0.000) 0.888 


yii 


(-0.193) -0.343 


(-0.193) -0.371 


(-0.193) -0.414 


yi2 


(1.158) 2.144 


(1.158) 2.296 


(1.158) 2.482 



Table 6: NLO Wilson coefficients at = 1 GeV in the HV scheme for = 170 GeV 
(a = 1/128). The corresponding values at = myj/ are given in parenthesis. In 
addition, at // = nic we have 2:3-12 (mc) = 0. The coefficients y\\ and y\2 are given at 
the "LO" (QCD two-loop mixing) and their values are 75-scheme independent. 
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e'/e X 10"* (NLO) 


rrif 


150 GcV 


170 GcV 


190 GcV 




1 (juad. 


11 (iuad. 


1 quad. 


11 (juad. 


1 quad. 


11 quad. 


Qs 


0.3 


4.1 


0.2 


2.9 


0.3 


3.8 


0.2 


2.6 


0.3 


3.5 


0.2 


2.3 


Qa 


-4.7 


-3.3 


-4.4 


-3.0 


-4.0 


-2.7 


Qs 


-0.6 


-0.4 


-0.5 


-0.4 


-0.5 


-0.3 


Qe 


9.1 


6.3 


8.4 


5.6 


7.8 


5.1 


Qt 


0.7 


-1.8 


0.5 


-1.2 


0.4 


-2.7 


0.2 


-1.8 


0.0 


-3.6 


0.0 


-2.4 


Qs 


-4.7 


-3.3 


-5.3 


-3.6 


-6.0 


-3.9 


Q9 


2.9 


2.0 


3.1 


2.1 


3.2 


2.1 


Qio 


-0.7 


-().."•) 


-0.8 


-0.5 


-0.8 


-O..") 


Qii 


-0.4 


-0.3 


-0.3 


-0.2 


-0.3 


-0.2 


All 


2.0 


1.4 


0.8 


0.5 


-0.4 


-0.2 



Table 7: Anatomy of e' je for Ag^^ = 200 MeV {as{mz)-j^ = 0.109), B^ ' > = 

(3/2') (1/2) 

B^' ^ = B^-j^ = 1, in the HV scheme. The contribution of each operator is shown at 
/i = 1 GeV, together with partial grouping of the gluonic and the electroweak sectors. 
The contribution of Q12 is being neglected. The entries for — Qio represent the 1 /Nc 
central values (for a discussion on the error bars sec rcf. [2], whose input parameters 
have been thoroughly assumed). The Qn amplitude is computed in the xQM, using 
the "LO" Wilson coefficient and setting {qq)G = —(275 MeV)^ as a typical value (see 
Table 1). 
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e'/e X 10* (NLO) 


mt 


150 GcV 


170 GeV 


190 GeV 


5km 


1 quad. 


11 quad. 


I quad. 


II quad. 


I quad. 


II quad. 


Qz 


0.5 


6.2 


0.3 


4.3 


0.5 


5.7 


0.3 


3.8 


0.4 


5.3 


0.3 


3.5 




-5.0 


-3.4 


-4.6 


-3.1 


-4.2 


-2.8 


gs 


-0.6 


-0.5 


-0.6 


-0.4 


-0.5 


-0.4 




11.3 


7.8 


10.4 


7.0 


9.6 


6.4 


Q-r 


0.6 


-3.2 


0.4 


-2.3 


0.3 


-4.3 


0.2 


-2.9 


0.0 


-5.3 


0.0 


-3.5 


Qs 


-5.9 


-4.1 


-6.7 


-4.5 


-7.5 


-5.0 


Q9 


3.0 


2.1 


3.1 


2.1 


3.3 


2.2 


Qio 


-0.9 


-0.6 


-1.0 


-0.7 


-1.0 


-0.7 


Qu 


-0.4 


-0.3 


-0.4 


-0.2 


-0.3 


-0.2 


All 


2.6 


1.8 


1.1 


0.7 


-0.4 


-0.2 



Table 8: Same as in Table 5 for Ag^^ = 300 MeV {as{mz)j^ = 0.116). 
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e'/e X 10* (NLO) 


mt 


150 GcV 


170 GeV 


190 GeV 


5km 


1 quad. 


II quad. 


I quad. 


II quad. 


I quad. 


II quad. 


Qz 


0.7 


8.7 


0.5 


6.0 


0.7 


8.0 


0.5 


5.4 


0.6 


7.4 


0.4 


4.9 




-5.2 


-3.6 


-4.8 


-3.2 


-4.4 


-2.9 


gs 


-0.7 


-0.5 


-0.7 


-0.4 


-0.6 


-0.4 




13.8 


9.6 


12.7 


8.6 


11.8 


7.8 


Q-r 


0.5 


-5.0 


0.4 


-3.4 


0.2 


-6.1 


0.2 


-4.1 


0.0 


-7.3 


0.0 


-4.8 


Qs 


-7.4 


-5.1 


-8.4 


-5.7 


-9.4 


-6.2 


Q9 


3.2 


2.2 


3.3 


2.3 


3.5 


2.3 


Qio 


-1.2 


-0.9 


-1.3 


-0.9 


-1.4 


-0.9 


Qu 


-0.4 


-0.3 


-0.4 


-0.3 


-0.4 


-0.3 


All 


3.3 


2.3 


1.5 


1.0 


-0.3 


-0.2 



Table 9: Same as in Table 5 for Ag^^ = 400 MeV {as{mz)j^ = 0.122). 
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Fig. 1 

Figure 1: The relevant diagrams for the tt"*" transition induced by Qn. The 

circled cross indicates the operator insertion. 




.-00-. 
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Fig. 3 



Figure 3: The relevant diagrams for the K 
cross indicates the Qq insertion. 



TT transition induced by Qq. The circled 






Fig. 4 



Figure 4: Building blocks for the hadronic matrix elements. The diagram on the top 
corresponds to (0| S75M \K~^(k)), whereas the two diagrams on the bottom contribute to 
{■K^ {p)\sd\K~^ (k)) . The circled cross indicates the quark current (or density) insertion. 



36 



References 

[1] G.D. Barr et al. (NA31 Coll.), Phys. Lett. B 317 (1993) 233; 
L.K. Gibbons et al. (E731 Coll.), Phys. Rev. Lett. 70 (1993) 1203. 

[2] A.J. Buras, M. Jamin and M.E. Lautenbacher, Nucl. Phys. B 408 (1993) 209. 

[3] M. Ciuchini, E. Franco, G. Martinelli and L. Reina, Nucl. Phys. B 415 (1994) 403; Phys. 
Lett. B 301 (1993) 263. 

[4] R.S. Chivukula, J.M. Flynn and H. Gcorgi, Phys. Lett,. B 171 (1986) 453; 
A.J. Buras and J.-M. Gerard, Nucl. Phys. B 264 (1986) 371; 
A.J. Buras, J.-M. Gerard and R. Riickl, Nucl. Phys. B 268 (1986) 16; 
W.A. Bardccn A.J. Buras and J.-M. Gerard, Phys. Lett. B 192 (1987) 138, 156; 
G. Buchalla, A.J. Buras and K. Harlandcr, Nucl. Phys. B 337 (1990) 313. 

[5] A. Pich and E. de Rafael, Nucl. Phys. B 358 (1991) 311. 

[6] S. Weinberg, Physica 96A (1979) 327; 

A. Manohar and H. Georgi, Nucl. Phys. B 234 (1984) 189; 
A. Manohar and G. Moore, Nucl. Phys. B 243 (1984) 55; 
J. Bijnens, H. Sonoda and M.B. Wise, Can. J. Phys. 64 (1986) 1; 
J. Bijnens, Nucl. Phys. B 367 (1991) 709. 

[7] D. Espriu, E. de Rafael and J. Taron, Nucl. Phys. B 345 (1990) 22. 

[8] J. Bijnens, C. Bruno and E. de Rafael, Nucl. Phys. B 390 (1992) 501. 

[9] S. Bertolini, M. Fabbrichesi and E. Gabrielli, Phys. Lett. B 327 (1994) 136. 

[10] A.J. Buras, M. Jamin, M.E. Lautenbacher and RH. Weisz, Nucl. Phys. B 370 (1992) 69, 

(Addendum) ibid. 375 (1992) 501; 

A.J. Buras, M. Jamin, M.E. Lautenbacher and RH. Weisz, Nucl. Phys. B 400 (1993) 37; 
A.J. Buras, M. Jamin and M.E. Lautenbacher, Nucl. Phys. B 400 (1993) 75. 

[11] M. Ciuchini, E. Franco, G. MartinelU, L. Reina and L. Silvestrini, Phys. Lett. B 316 

(1993) 127; 

M. Ciuchini, E. Franco, L. Reina and L. Silvestrini, Nucl. Phys. B 421 (1994) 41. 

[12] A.J. Buras, M. Misiak, M. Munz and S. Pokorski, Munich preprint MPI-Ph/93-77; 

M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. Lett. B 334 

(1994) 137. 

[13] G. Ecker, H. Neufeld and A. Pich , Phys. Lett. B 278 (1992) 337; preprint CERN- 
TH.6920/93. 

[14] J. Bijnens, G. Ecker and A. Pich, Phys. Lett. B 286 (1992) 341. 

[15] J. Kambor, J. Missimer and D. Wyler, Nucl. Phys. B 346 (1990) 17; 
G. Ecker, J. Kambor and D. Wyler, Nucl. Phys. B 394 (1993) 437 

[16] J. Bijnens, Int. J. Mod. Phys. A 8 (1993) 3045. 



37 



[17] R. Alkofer and H. Reinhardt, preprint UNITU-THEP-19/1992 (unpublished); 
R. Ball and G. Ripka, preprint CERN-TH.7122/93; 
J. Bijnens and J. Prades, Phys. Lett. B 320 (1994) 130. 

[18] J. O. Eeg and I. Picek, Phys. Lett. B 301 (1993) 423. 

[19] J. O. Ecg and 1. Picek, Phys. Lett. B 323 (1994) 193. 

[20] C. Bruno and J. Prades, Z. Physik C 57 (1993) 585. 

[21] P. Damgaard and R. SoUacher, Phys. Lett. B 322 (1994) 131. 

[22] J.F. Donoghue and B.R. Holstein, Phys. Rev. D 32 (1985) 1152. 

[23] G. Feinberg, P. Kabir and S. Weinberg, Phys. Rev. Lett. 3 (1959) 527. 

[24] C. Bernard et al., Phys. Rev. D 32 (1985) 2343. 

[25] A. Di Giacomo, H. Panagopoulos and E. Vicari, Nucl. Phys. B 338 (1990) 294. 

[26] J.A. Cronin, Phys. Rev. 161 (1967) 1483. 

[27] H.-Y. Cheng, Int. J. Mod. Phys. A 4 (1989) 595. 

[28] J.F. Donoghue, Phys. Rev. D 30 (1984) 1499 

[29] Y. Dupont and T.N. Phani, Phys. Rev. D 29 (1984) 1368; 
T.N. Pham, Phys. Lett. B 145 (1984) 113. 

[30] M.B. Gavela, A. Le Yaouanc, L. Oliver, O. Pene and J.C. Raynal, Phys. Lett. B 148 

(1984) 253. 

[31] J. Gasser and H. Leutwyler, Ann. Phys. (NY) 158 (1984) 142. 
[32] S. Fajfer, preprint CERN-TH.7096/93. 

[33] S. Bertolini, F. Borzumati and A. Masiero, Phys. Rev. Lett. 59 (1987) 180; 

N.G. Deshpande, P. Lo, J. Trampetic, G. Eilam and P. Singer, Phys. Rev. Lett. 59 (1987) 
183; 

B. Grinstcin, R. Springer and M. Wise, Phys. Lett. B 202 (1988) 138; Nucl. Phys. B 
339 (1990) 269; 

R. Grigjanis, P.J. O'Donnell, M. Sutherland and H. Navelet, Phys. Lett. B 213 (1988) 
355, (E) ibid. 286 (1992) 413; 

G. Cella, G. Curci, G. Ricciardi and A. Vicere, Phys. Lett. B 248 (1990) 181; 
M. Misiak, Phys. Lett. B 269 (1991) 161; Nucl. Phys. B 393 (1993) 23; 
K. Adcl and Y.P. Yao, Mod. Phys. Lett. A 8 (1993) 1679; 
M. Ciuchini et al., as in ref [11]. 

[34] L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rep. 127 (1985) 1. 



38 



